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A possible definition of strong/symmetric hyperbolicity for a second-order system of evolution 
equations is that it admits a reduction to first order which is strongly /symmetric hyperbolic. We 
investigate the general system that admits a reduction to first order and give necessary and sufficient 
criteria for strong/symmetric hyperbolicity of the reduction in terms of the principal part of the 
original second-order system. An alternative definition of strong hyperbolicity is based on the 
existence of a complete set of characteristic variables, and an alternative definition of symmetric 
hyperbolicity is based on the existence of a conserved (up to lower order terms) energy. Both these 
definitions are made without any explicit reduction. Finally, strong hyperbolicity can be defined 
through a pseudo-differential reduction to first order. We prove that both definitions of symmetric 
hyperbolicity are equivalent and that all three definitions of strong hyperbolicity are equivalent (in 
three space dimensions). We show how to impose maximally dissipative boundary conditions on any 
symmetric hyperbolic second order system. We prove that if the second-order system is strongly 
hyperbolic, any closed constraint evolution system associated with it is also strongly hyperbolic, 
and that the characteristic variables of the constraint system are derivatives of a subset of the 
characteristic variables of the main system, with the same speeds. 
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I. INTRODUCTION 

Research in numerical relativity has recently focused 
on obtaining a well-posed continuum initial-boundary 
value problem as a starting point for numerical time evo- 
lutions of systems such as a black-hole binary. Well- 
posedness of an initial-boundary value problem implies 
that an estimate 
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VI. The evolution of constraints on the original second-order svstemfe 

exists, where u(x,t) is the solution 
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VII. Symmetric hyperbolic-parabolic svstemd 10 



VIII. Conclusiond 



u(x, 0) the initial 

data, g(x,t) appropriate free boundary data, S denotes 
a linear perturbation and || • || stands for appropriate 
norms (which may involve spatial derivatives), and where 
F(t) is independent of the initial and boundary data. 
This means that the solution depends continuously on the 
initial and boundary data. Hyperbolicity is a property of 
t he evolution equations that can be used as an algebraic 
Diagonalisabilitv of matrices with a zero blockt tiliterion for well-posedness. We briefly review several 

notions of hyperbolicity. 
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Consider a system of quasilinear evolution equations 
that is first order in both space and time, or 

u = P l (u)u^ + S(u), (2) 

where u is a vector of variables and P % are square matri- 
ces. 

Definition 1: The system 0) is called weakly hyperbolic 
if the matrix P n = riiP % has real eigenvalues for any unit 
vector ni . 

Definition 2: The system 0) is called strongly hyper- 
bolic if P n is diagonalisable with real eigenvalues for any 
ni, and the matrix T n that diagonalises it and its inverse 
T^ 1 depend smoothly on n^. 

Definition 3: The system 0) is called symmetric hyper- 
bolic if there exists a Hermitian, positive definite matrix 
H such that HP n is Hermitian for any direction rii and 
where H does not depend on Hi . 

The following properties of strongly and symmetric hy- 
perbolic systems give a more practical meaning to the 
definitions, and we shall use them later to define strong 
and symmetric hyperbolicity for second-order systems. 
The key concept for strong hyperbolicity is 
Definition 4: A characteristic variable with speed —A in 
the rii direction is a linear combination u of the variables 
u that obeys 

d t u = Xd n u + (3) 

where rii is normalised with respect to some metric, d n = 
n l di, and the dots denote derivatives transverse to rii with 
respect to the same metric, and lower order terms. 

If we write uasu = u*u, where u is a constant vector 
of coefficients, then 

d t (v)u) = v)P n d n u + . .. = \d n (v)u) + .. . (4) 

if and only if Tv is a left eigenvector of P n or equivalently 
if u is an eigenvector of P* 1 '. Characteristic variables u 
of the first-order reduction therefore correspond to left 
eigenvectors of P n . This gives us 
Lemma 1: A first- order system is strongly hyperbolic 
if and only if it admits a complete set of characteristic 
variables with real speeds that depend smoothly on rii. 

The key concept for symmetric hyperbolicity is that of 
an energy: 

Definition 5: An energy e is a quadratic form in u that 
is positive definite in the sense that e = if and only 
if u = 0, and which is conserved in the sense that there 
exists a flux <j) 1 quadratic in u such that 

e = 0V (5) 

With 

e = vtHu, ft = u f HP l u, (6) 

we have 

Lemma 2: A linear first- order system with constant co- 
efficients is symmetric hyperbolic if and only if it admits 
an energy. 



For quasilinear systems this energy is conserved in the 
approximation where S(u) is neglected and P l (u) is ap- 
proximated as constant. (Physically, this corresponds 
to considering small high-frequency perturbations 5u.) 
When boundaries are present, the time derivative of the 
energy can be estimated in terms of free boundary data. 

Strong hyperbolicity of a first-order system is nec- 
essary and sufficient for a well-posed Cauchy problem. 
The Cauchy problem for a merely weakly hyperbolic sys- 
tem is typically ill-posed in the presence of lower-order 
terms. Symmetric hyperbolicity implies strong hyperbol- 
icity, and is therefore also sufficient for well-posedness of 
the Cauchy problem. Furthermore, symmetric hyperbol- 
icity can be used to prove well-posedness of the initial- 
boundary value problem for a certain class of boundary 
conditions called maximally dissipative [lj. 

Hyperbolicity for equations or systems of equations of 
higher than first order is less well-established. A defini- 
tion of weak hyperbolicity exists for systems of arbitrary 
order, but as for first-order systems, it does not guaran- 
tee well-posedness Q • Alternatively, a quasilinear system 
that is second order in both space and time, or 

P^(u, du)u„ v + S(du, u) = 0, (7) 

is called hyperbolic if P^ is a Lorentzian metric, that is 
if the principal part of the system is that of a wave equa- 
tion [3( ■ Christodoulou has recently generalized the idea 
introducing the concept of regular hyperbolicity, with less 
strict positivity requirements on the elliptic block of the 
principal part 0. Both can be used as criteria for well- 
posedness of the Cauchy problem. The Einstein equa- 
tions are second order, but they fit these definitions of 
hyperbolicity only when written in harmonic gauge, and 
so there are no standard definitions of hyperbolicity im- 
mediately applicable to forms of the Einstein equations 
commonly used in numerical relativity. 

One possible approach to the well-posedness of a 
second-order system is to reduce it to first order by intro- 
ducing auxiliary variables, and to define the second-order 
system to be strongly hyperbolic or symmetric hyperbolic 
if the reduction is. An ad-hoc definition alongthose lines 
has been used by Sarbach and co-authors [oj. |G| to prove 
well-posedness of the BSSN formulation of the Einstein 
equations. We formalise this approach in Section ITm 

Independently, Nagy, Ortiz and Reula 7], following 
Kreiss and Ortiz y| have used a pseudo-spectral reduc- 
tion to define strong hyperbolicity. This method does 
not appear to generalise to symmetric hyperbolicity, in- 
tuitively because Fourier transforms cannot be carried 
out on a domain with arbitrary boundary. We briefly 
review this approach in Section By casting it in the 
notation of Section ITTll show that the two definitions of 
strong hyperbolicity are equivalent. 

As a third alternative, Gundlach and Martfn-Garcfa 
P.IIo| define strong and symmetric hyperbolicity directly 
from the second-order system, by focusing on the exis- 
tence of characteristic variables in strong hyperbolicity, 
and of an energy in symmetric hyperbolicity. We review 
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this approach in Section llVl and show that its definitions 
of both strong and symmetric hyperbolicity are equiva- 
lent to those using a first-order reduction. 

Outside the main line of this paper, we analyse in Sec- 
tion lVII the well-posedness of the propagation of any con- 
straints that the original second-order system is subject 
to. In Section IVIII we apply our results for symmet- 
ric hyperbolic systems to mixed symmetric hyperbolic- 
parabolic systems. Section IVIIII summarises our results. 



II. THE SYSTEM 

In this short Section, we establish notation for the class 
of system that we want to investigate, and clarify the 
relation between systems that are second order in both 
time and space, only in space, or only in space and that 
only in some of the variables. 

We begin systems that are first order in time. Formula- 
tions of the Einstein equations based on the ADM formu- 
lation are naturally first order in time and second order 
in space. Some cannot even be written in second-order in 
time form (for example because they have an odd num- 
ber of variables). With a non-zero shift, the first-order 
form may also be preferable for numerical simulations 
|llj |. For simplicity, we restrict attention to linear sys- 
tems with constant coefficients. These can be considered 
as the linearisation and frozen coefficients approximation 
of a nonlinear system. 

The class of systems that we are interested in are not 
uniformly second-order: some variables u may appear in 
the evolution equations without second spatial deriva- 
tives. Writing u — (v,w), where w are those variables 
that appear only with first derivatives, we consider there- 
fore 

v = A^v^j + A\v ti + Aiv + A 2 w ti + A 2 w + a, (8) 
ib = Bfv ti j + B{v,i + B x v + B\w,i + B 2 w + b. (9) 

Here v and w are column vectors (not necessarily of the 
same length) of variables, the capital letters represent 
constant matrices of the appropriate dimension, and a, 
b are forcing functions. We assume that while at least 
one second derivative of every variable v appears in the 
equations, the number of variables v has been minimised. 

To reduce the system <|8I9|) to first order, we define the 
auxiliary variables di = Vj . By taking a spatial derivative 
of JSJ|, we find 

di = A{ k v tij k+A{v,ij +A 1 v ti + A 3 2 w t ij + A 2 w ti + a ti . (10) 

This is of a higher order than we started from, unless 
AY and A 2 both vanish. (An example of a second-order 
evolution equation that cannot be reduced to first order 
is the heat equation u = u" .) We have 
Lemma 3: The general second-order in space, first-order 
in time linear system that can be reduced to first order by 



the introduction of auxiliary variables is of the form 

v = A\v t j + Aiv + A 2 w + a, (11) 

w = />•; •'<•., ,, + B\v. % + Biv + B l 2 w^ + B 2 w + b. (12) 

From now on we refer to this as "the" second-order sys- 
tem. We have underlined the highest derivatives. With- 
out loss of generality we assume from now on that B^ is 
symmetric. 

In order to understand how general the system (|11I12|) 
is, it is interesting to convert it into second-order in both 
space and time form. Taking a time derivative of l|ll|) 
and using (|12|) to replace ib, we obtain 

v = A 2 B[ J v, tj + A\v ti + A 2 B l 2 w A + A 2 B 2 w + ... (13) 

where we have written out all second derivatives and all 
appearances of w. We can eliminate the remaining ap- 
pearances of w and w t i in terms of v using i|ll|) if and 
only if the matrices of the system obey 

rank(A 2 ) = rank ( ^ ) = rank ( ^ ) . (14) 

When A 2 is invertible, which in particular implies equal 
numbers of v and w variables, these conditions are auto- 
matically obeyed. On the other hand, any fully second- 
order system in a set of variables v can be reduced to 
the form (|11I12|) by introducing v = w. Therefore the 
class of first-order in time, second-order in space systems 
(|11I12|) includes the class of fully second-order systems, 
but is much bigger. 

III. FIRST-ORDER REDUCTION METHOD 

A. Parameterised reduction 

In reducing i|llll2fl to first order by defining di = v t i, 
we can write each occurrence of v t i also as di, or a mix- 
ture of the two, and similarly we can write v^j as dij or 
djj. To parameterise these ambiguities, we formally add 
multiples of the auxiliary constraint 

c 4 ee di - v t i = (15) 

and its antisymmctrised derivative 

°ij = d U-i] = C [jA = ( 16 ) 

to all three equations. We could not add cuj\, or any 
higher derivatives of the auxiliary constraints, without 
increasing the order of the system. The general reduction 
to first order is therefore 

v = A\v t i + A\v + A 2 w + a 

+Ai j c ij +A i 3 c i , (17) 
lb = B\ 3 dij + B{v,i + Biv + B l 2 w. t + B 2 w + b 

+BYc ij + Blc i , (18) 
di = A[dij + A\v.i + A 2 w.i + a t i 

• n!c h . ■ !)/"<■,,. (19) 
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From now on, we shall refer to this system as "the" reduc- 
tion. Wc shall refer to the constant matrices A l 3 , A 3 , B\, 
B 3 , IV and ZV as the reduction parameters. Without 
loss of generality, we assume that A 3 , B 3 and Dk J are 
antisymmetric in i and j. The terms in the second-order 
system that become principal terms in the reduction are 
the ones underlined in l|lH12fl . Wc shall call these the 
principal part of the second-order system. 



Definition 3a: The second-order system is defined to be 
symmetric hyperbolic if and only if it admits at least one 
reduction to first order that is symmetric hyperbolic. 

In the remainder of this Section we derive necessary 
and sufficient conditions for these definitions to hold, for- 
mulated directly in terms of the principal part of the 
second-order system, without reference to a reduction. 



B. Auxiliary constraint evolution 

The evolution of v and implies an evolution of the 
auxiliary constraints Cj. This auxiliary constraint sys- 
tem can be written in first-order in space and time form 
by introducing Cjj = cua (as already defined above) as 
auxiliary variables. This results in 



A\ci, 



A 3 Cj t i 



A lk r , 
^3 C J*V 



.) k ■ 



(20) 
(21) 



As the right-hand side of this system of linear PDEs is 
homogeneous, c(x,0) = implies c(x 1 0) = 0. Assuming 
that the coefficients are constant, as we do in this paper, 
on taking a Fourier transform in x l we obtain a sepa- 
rate ODE for each wavenumber uj 1 , and it follows that 
c(x, t) = is the unique solution with c(x, 0) = 0. 

We have shown that if the auxiliary constraints are 
zero initially, they remain zero at all times. This allows 
us to prove well-posedness for the reduction to first order, 
and then restrict it to the subset of solutions that obey 
the auxiliary constraints in order to infer well-posedness 
of the original second-order system. 

Note that in order to make this argument, well- 
posedness of the auxiliary constraint system is not re- 
quired, as we require only existence and uniqueness of 
the zero solution, not estimates of any non-zero solution. 



D. 2+1 split 

We introduce a matrix notation in the groups of vari- 
ables v, w and d. In this notation, ll'L'li is 





(23) 



where 



P* 



A\-A\ 
B\-B\ 
AiS/ - Dj 





B% 



Bf 



Aik 
A 3 



AW 



Bf 
-IV 



(24) 



We expand the tensor indices i, j and k in 1)24(1 into 
their components in the direction ni and transverse to 
it. For this purpose we need a positive definite metric 
tensor 7.^ . This can be chosen to be Sij , or the physical 
3-metric for example in applications to general relativity. 
We can then define the normal component of a tensor 
such as d n = n l di where n % = 7 y 7ij and ni is normalised 
so that niUjj v = 1, and its transverse part such as (5^ — 
nin J )dj, which we denote by dA- For a reason that will 
become apparent immediately, we also re-order the rows 
and columns. We call the resulting matrix V . It is related 
to P n by a unitary transformation that depends on n^, 
but it is not the component of any vector of matrices. 
We have 



C. Definition of hyperbolicity 

We now focus on the principal part of the first-order 
reduction, 



d t u ~ P l diU, 



(22) 



where here and in the following ~ denotes equality up to 
lower-order terms and now u stands for (v,w,di). The 
degree of hyperbolicity of the first-order system depends 
on the the reduction parameters. The appropriate defi- 
nitions of hyperbolicity are therefore the following: 
Definition la: The second-order system is de- 

fined to be weakly hyperbolic if and only if it admits at 
least one reduction to first order \11%19\) that is weakly 
hyperbolic. 

Definition 2a: The second-order system is defined to be 
strongly hyperbolic if and only if it admits at least one 
reduction to first order that is strongly hyperbolic. 



/ w \ 




/ w 


d n 




d n 






V 


\d A ) 




\d B 



where 



V 



( B% 
A 2 


V 



A\ 





Bl 



Bf 
Ax 
A'l 
-D A n 



transverse derivatives, 

(25) 



3 

LJn 
An 
^3 



B nB 

D nB 

Un 
A 3 



A?S A * + D A nB ) 



(26) 



We write this in shorthand form as 



V = 



A B 
C 



A = 



B% 
A 2 



A\ 



(27) 



The eigenvalues of V are those of P n , and one is diag- 
onalisable if and only if the other is. Therefore we can 
investigate weak and strong hyperbolicity using V . The 
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fact that V has a zero block for all choices of the reduc- 
tion parameters has several interesting consequences for 
its eigenvalues and eigenvectors, which respectively de- 
termine the weak and strong hyperbolicity properties of 
the first-order system. 



E. Eigenvalues of V: weak hyperbolicity 

The first of these consequences is that the eigenvalues 
of V (and hence of P n ) are the union of the eigenvalues 
of A and C, and are independent of B. We will show in 
Section lTTTFl that we can choose the reduction parameters 
so that C has real eigenvalues. This gives us 
Lemma 4: The second-order system is weakly hyper- 
bolic according to Definition la if and only if A has real 
eigenvalues for all rii . 

We can further analyse the eigenvalues of A. If we re- 
place rii by —ni, then B^ and A r { change signs, while B\ n 
and Ai do not. This means that if A is an eigenvalue of A 
for some rii then — A is an eigenvalue of A for — n^. The 
eigenvalues must therefore either be A = O(rii) where O 
is odd in m, or they are paired as X± = O(rij) ± E(rii) 
where O is odd and E is even. As a consequence, if the 
dimension of A is odd, at least one eigenvalue must be 
of the form A = 0(rii), and by continuity it must vanish 
for some m. 



F. Eigenvectors of V: strong hyperbolicity 

The reduction is strongly hyperbolic if and only if P n 
(or equivalently V) is diagonalisable. In Appendix 1X1 we 
show that a necessary condition for V to be diagonalis- 
able is that both A and C are diagonalisable. If the sets 
of eigenvalues of A and C are disjoint, this is also suffi- 
cient. If they have any eigenvalues in common, additional 
necessary criteria arise which involve B. 

To state these conditions in a simple form, diagonalise 
A and C simultaneously, so that 



S~ l VS 



A.4 B 
A c 



(28) 



where A ^4 and Ac are diagonal matrices with the eigen- 
values of A and C respectively in the diagonal. Let all 
repeated eigenvalues be grouped together in these ma- 
trices. Then, for each eigenvalue common to A and C, 
the corresponding block of B must vanish if V is to be 
diagonalisable. 

The block A does not contain any reduction parame- 
ters, and so is determined by the original second-order 
system. For a given direction rii, blocks B and C are de- 
termined completely by the choice of reduction parame- 
ters, but there are not enough reduction parameters to 
make this true for all directions ni at once. (For exam- 
ple, as B-y is symmetric and B§ is antisymmetric in ij, 
B'i B + B$ B can be made to vanish for any one direction 
rii, but not for all directions.) 



We shall consider one choice of reduction parameters 
in discussing strong hyperbolicity (here) and another one 
in discussing symmetric hyperbolicity (in the next sub- 
section). Both choices set 



A 10 

4 



= 0, 

= A\, 

= B\, 

= AiSi 3 , 



(29) 



This partial choice has the effect of decoupling v from 
the w and di. To discuss strong hyperbolicity in the 
case of three space dimensions, we complete the choice 
of reduction parameters by 



-°3 



D 



i k 



o, 

5j i A\ - 6j k A\ + ifief 



(30) 



where fi is a real constant, e^-fe is the totally antisym- 
metric tensor in three dimensions and i = This 
gives 



I B 



V 



\ 



j B[ ln B[ lB 

A 2 A1 Af 



i\xt A nB ) 



(31) 



The diagonal block ifieA nB is diagonalisable with real 
eigenvalues ±/i. The complexification is unusual, but 
the term multiplied by i[i is proportional to the auxil- 
iary constraint Cy , and so has no influence on the orig- 
inal second-order system. If we choose \x large enough, 
the eigenvalues of this block are distinct from those of 
the complementary diagonal block (containing A and a 
zero row and column), and we have found a choice of re- 
duction parameters that makes V diagonalisable if A is 
diagonalisable. The existence of this choice of reduction 
parameters completes our proof of 

Theorem 1: The second-order system is strongly hy- 
perbolic according to Definition 2a if and only if A is 
diagonalisable for all rii where the diagonalising matrix 
and its inverse depend smoothly on rii. 

Our proof of this theorem assumes three space dimen- 
sions, but we suspect that the theorem holds in any num- 
ber of space dimensions. 

Note that this criterion is based only on the coefficients 
on the second-order system, without explicit reference to 
the reduction. Note also that for the choice of reduction 
parameters we have used here, the auxiliary constraint 
system is strongly hyperbolic, see Appendix 151 



G. Symmetric hyperbolicity 

According to Definition 3, the reduction is symmetric 
hyperbolic if and only if there is a Hermitian matrix H 
such that 



(iJP") t = HP" 



(32) 
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for all rii, with H independent of rii and positive definite. 
Note that in the definition (|32[l we cannot replace P n by 
V. 

We again make the partial choice (|29[l of reduction 
parameters. B^ and £V will be determined in the fol- 
lowing in terms of H. The resulting form of P % is [see 
(|2"3j) for the definition of P l \ 



P' 






B\ k + B\ 



A>d/ A\Sj + Dj l 



(33) 



Clearly it is sufficient to find a symmetriser for the (w, di) 
block. We parameterise H as 



(34) 



with K and M Hermitian and positive, IO = K, K > 0, 
and M* mj = M jm , M > 0. 

The nontrivial, (w,di) block of HP 1 is 



f 












K 




w 







M m 3 J 





f KB\ + UA 2 
lt m B\ + M mi A 2 

V 



K(B\ 



ik 



Bf) 



\ 



+L k A\ +LW/ 
-M mk A[ + M m W/ k ) 



(35) 



A necessary condition for this to be Hermitian is for the 
matrix 



KB% + L n A 2 KB™ n + L n A^ 
tf n BH + M nn A 2 L^ n B" n + M nn A n 



(36) 



to be Hermitian for all rii. This is just the condition that 
A admits a symmetriser, or 



HA={HA)\ where H 



jj[n ]\/[nn ) i (37) 



for all rii. If the system is strongly hyperbolic, A always 
admits a symmetriser formed from its eigenvectors. Nev- 
ertheless, (|37|) is a non-trivial condition because Tt must 
form a part of H, so that its blocks L n and M nn are 
given by L l rii and M^riiUj, where K, L l and M y do 
not depend on rii. We now show that (|37|) actually im- 
plies that all of (|35|l is Hermitian for a particular choice 
of reduction parameters B$ and Dj . We do this for 
each block of (|35|) in turn. 

From the top left block of l|35|l we have 



(38) 



[KB\ + UA 2 y = KB l 2 + L l A 2 



This equation does not contain any reduction param- 
eters, and it is clearly equivalent to the condition of 
KB 2 + L n A 2 being Hermitian for all rii, which is con- 
tained in the first diagonal block of Ij37(l . 
The off-diagonal blocks of (|35fl give 

[L^ k B l 2 + M ki A 2 ) ] = K(B\ k + Bf) + L k A\ + L j D 3 lk . 

(39) 



If we denote this equation by T lk = 0, then its symmetric 
part = 0, or 



B\ [l L k) + A\h& k) = KB\ k , +L {k A\ 



(40) 



again does not contain any reduction parameters. Fur- 
thermore = if and only if T ik n t n k = for all 
rii (this is a property of all totally symmetric tensors), 
and this is precisely the condition contained in the off- 
diagonal terms of l|37|) . The antisymmetric part = 0, 
or 



Bt [i L H l + A\m^ = KB lk + L j Dj tk , 



(41) 



can be solved for B$ because K is by assumption invert- 
ible. 

Finally, the bottom right block of (|35|l gives 



L imr B fk + B f) + M mk A\ + M mj D, 



ik 



= [tf k {Bf 1 + B l 3 m ) + M km A[ + M k W/ m ] f . (42) 

If we write this as T mlk = 0, the totally symmetric part 
T (mik) = 0; Qr 

L Km B ik) + M (mfe A i) = \ L Kk B ^) + M ( km A]>} f . (43) 

does not contain any reduction parameters. It is equiv- 
alent to T mlk n m riink = for all rii, and so vanishes 
because of the last diagonal block of i|37|) . After solving 
(gTJ for Bz j we write (g2J) as 

jjmik + Mmj D .ik = jjikim + ( \/ /■.-/;/»• ) j (44) 

where we have defined 

jjmik = i) m B f + M mk A\ 

+L im K -i f B ili L k] + AlM^ ik A , (45) 
M ij = M ij - L^K- X V. (46) 

M is Hermitian (M™ = M 3% ) and positive definite [be- 
cause H is positive definite in particular when restricted 
to vectors (v,w,di) = (0, — K~ 1 L h dk, di)], and hence in- 
vertible. We define 



j^mik _ jj\kim jjmik 



D mik = M mj D .ik^ 



and so 

while T( mlk ) = is equivalent to 
(|4^|l becomes 

jymik jjjkim j^mik 



(48) 



(49) 



(50) 
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which has the general solution 



variable in the direction m because 



6D mik = X kml + 2X klm + 3X 



ikn 



_^j£imk _|_ ^J^ m *^ _|_ Y 771 '^ 



(51) 



d t (v,A) = d A (d t v) 



(53) 



using both JJSJ and (gHJ. The object Y mik must obey 

^mik -y^fmifc] yrnik y^rnik (52) 

but is otherwise arbitrary. It parameterises the part of 
the Dj lk that is not determined by H, and can be set to 
zero. 

We have shown that a necessary and sufficient condi- 
tion for P % to admit a symmetriser H is for A to admit a 
symmetriser Ti that depends on rii in the particular way 
given in 1|37|) . Therefore we have 

Theorem 2: A necessary and sufficient condition for the 
second-order system to by symmetric hyperbolic according 
to Definition 3a is that \3 r J\j holds for all rii for some 
H > parameterised by {34M - 

Note that positive definiteness of Ti does not imply 
that of H, because positive definiteness of M nn for all rii 
does not imply positive definiteness of M. The difference 
can be expressed in standard terminology as follows: a 
double quadratic form M lJ ab is rank-1 positive if and 



only if M %3 ABniUjmAmP 



and 



^nf > for all n 



it is 

A 



> for all n 
rank-2 positive if and only if ABn A 
Rank-2 positivity implies rank-1 positivity, but they are 
not equivalent when the indices A, B belong to spaces 
of dimension 3 or larger, as shown in the example given 
in 0| . This suggests it may be useful to introduce an 
intermediate concept of hyperbolicity based on positivity 
of Ti for all rii, rather than positivity of H. This already 
guarantees that A is diagonalisable, so that the system is 
strongly hyperbolic. The imposition of rank-1 positivity 
is also the key ingredient of the definition of "regular 
hyperbolicity" by Christodoulou Q for second order in 
both space and time systems, which is also known to yield 
well-posed problems. We believe there is a connection 
between regular hyperbolicity and our condition Ti > 0, 
but have not been able to show it. 



IV. DIRECT SECOND-ORDER METHOD 
A. Strong hyperbolicity 

Following [Tol | , we now elevate Lemma 1 to a definition 
of strong hyperbolicity for second-order in space, first- 
order in time systems, with the only difference that the 
u become linear combinations of v i and w: 
Definition 2b: The system is called strongly 

hyperbolic if and only if is there is a set of characteristic 
variables u linear in w and v t i obeying with real speeds 
A and where the matrix relating u and u and its inverse 
depend smoothly on ni . 

This definition has two interesting consequences. The 
first is that v a can always be considered as a zero speed 



so that the right-hand side can be considered as a sum 
of transverse derivatives only. (In fact arbitrary linear 
combinations of the v : a can be given arbitrary speeds.) 

The second consequence is that nontrivial character- 
istic variables are unique only up to adding transverse 
derivatives. If u is a characteristic variable with speed 
— A then, for any vector X A (made from the u) 



d t (u + d A x A ) 
d t u + d A (d t x A ) 

Xd n u + transv. deriv. 

\d n u + \d A (d n x A ) + 



transv. deriv. 



A<9 rl (u + 8aX a ) + transv. deriv. 



(54) 



and so u + 8aX a is also a characteristic variable with 
speed —A. Such calculations rely on commuting partial 
derivatives to interpret dAd n . . . either as a normal or a 
transverse derivative, depending on the situation. 

The second-order system has no reduction parameters. 
However, for the purpose of comparing the second-order 
approach with the reduction approach, we can translate 
the second-order approach into the notation of a first- 
order reduction. We account for the fact that and 
di, and v^j and vji, are now identical, by allowing the 
reduction parameters to depend explicitly on rii (so that 
they are not tensors.) This allows us to set the blocks B 
and C of V to arbitrary values for every n, , and this allows 
us to make arbitrary linear combinations of v t A and v 
characteristic variables with arbitrary speeds, and to add 
arbitrary combinations of v. a and v to any characteristic 
variables made from w and w n , as we have discussed 
above. A being diagonalisable is a necessary condition 
for V to be diagonalisable, and with B = and C = 
it is also sufficient. Therefore Definition 2b is equivalent 
to A being diagonalisable and by Theorem 1 it is then 
equivalent to Definition 2a. 

Alternatively, we can write down the principal part of 
the second-order system in a 2+1 split as follows: 

v ~ 0, (55) 
»,» - A[(v tn ) tj + A 2 w, n , (56) 
v.a ^ A[(v t j) tA + A 2 w <A , (57) 

+B ab {v, a ),b+B 1 2 w, u (58) 

where v i on the right-hand side is now considered a lower 
order term. Note the way the second derivatives have 
been ordered differently in u „ and v t A- In the language 
of reduction this corresponds to the reduction parameters 
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depending explicitly on rii. The matrix V becomes 



/ w \ 




( w \ 










+ transverse derivatives 










\v, A J 









where 



V 



B$ Bf n 

A 2 A™ 







(59) 



(60) 



Once again V is diagonalisable if and only if A is. Either 
way, and taking into account the smoothness conditions, 
we have 

Theorem 3: Definition 2b is equivalent to Definition 2a. 



B. Symmetric hyperbolicity 

Following , we now elevate Lemma 2 to a definition 
of symmetric hyperbolicity for second-order systems in 
Definition 3b: The second-order system is called sym- 
metric hyperbolic if and only if it admits a positive defi- 
nite energy e and a flux (j) 1 , both quadratic in w and Vj, 
that obey 

e-^V (61) 
Note that in obtaining (|61[) one can make arbitrary use 

Of Ujj = Vji. 

We parameterise the energy e by l|34|) , with di replaced 
by v t i, and parameterise the flux </> l as 



0* = (i«t,t>t ) 



F 



pik 

<mik 



W 



(62) 



with = F l and F^ klm = F mik . In the second-order 
system, there are no reduction parameters, and so the 
non-trivial part of HP 1 is given by 



KB\ + UA 2 



KB? + L k A\ 



L^ m B\ + M mi A 2 L^ m B\ k + M mk A\ 



(63) 



In a first-order system, energy conservation is precisely 
equivalent to HP 1 being Hermitian. In the second-order 
system, the relation between HP % and energy conserva- 
tion is more complicated, because = v t ki- 
To see this more clearly, we write out 



] {KB 2 + L l A 2 )w, 



-w\KB\ k + L^ k A^)v^ k 
+vl m (L^ m B t 2 + M m A 2 )w l 

+v} m {L^ m B\ k +M m ^ k A l })v, lk 

w^F l w tl +w ] F {lk) v tlk 
+v^ im Wi + v] m F m ^v, ik 



(64) 



(65) 



Comparing the first terms in e and (j) 1 t i we have 



F' 
F i 



KB\ + L l A 2 , 



(66) 
(67) 



which admits a solution F l if and only if l|38|) is obeyed. 
Comparing the second and third terms we have 



piik) = KBf+L^ k A^, 

pfun = L im B i +M mi A 2 , 



(68) 
(69) 



which admits a solution F lk if and only if (|40(l is obeyed. 
Comparing the fourth terms we have 

pm(ik) _ pfmj^ik _|_ j^,pn(k _ gmik (70) 
p\kim pmik (71) 

These admit a solution F mlk if and only if l|43|) is obeyed. 
It is clear that l|43[l is necessary. To demonstrate that it 
is sufficient we solve explicitly for F m%k . The general 
solution of (|70() is 

pmik gmik _j_ jmik (72) 

where J mik = J m i lk \ . Defining 

ymik _ £t~\kim gmik (73) 

the general solution of (|71fl is 

Qjmik ykmi _j_ c^ykim _j_ ^yikm 

+4: yimk + 5 ymik + yymik ^ (j^ 

if and only if yi mik ) = 0, which is equivalent to g3J|. The 
remaing free coefficient W mlk = W^ mlk ^ = W^ mik pa- 
rameterises terms in (j? whose divergence vanishes iden- 
tically. It can be set equal to zero without loss of gener- 
ality. 

We have shown that the second-order system admits a 
conserved energy if and only if i|38|) . I|40|) and (|43[) hold, 
which together are equivalent to (|37l) . This is equivalent 
to the existence of a conserved energy for the first-order 
reduction. The two energies are in fact the same un- 
der the (unambiguous) identification of di with v i. This 
means that the two definitions of symmetric hyperbolic- 
ity are equivalent, as both are equivalent to the matrix 
H defined by (|3~4*Jl being positive definite, and obeying 
(|37jl for all ni. We have 

Theorem 4: Definition 3b is equivalent to Definition 3a. 

Given an energy H of the second-order system, a first- 
order reduction that admits the same energy is given by 
the reduction parameters l|29|) and B^ and Dj deter- 
mined in Section IlII Gl Going the other way, an energy 
H for any first-order reduction is clearly also an energy 
for the original second-order system. If the second-order 
system admits a multi-parameter family of energies, then 
some of these parameters define reduction parameters 
and Dj of the first-order system, and the remain- 
der parameterise the energy of that particular reduction. 
An example of this split is given in Appendix iDl 
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C. Maximally dissipative boundary conditions 

A first-order symmetric hyperbolic system on a domain 
fl admits an energy 

E= I edV (75) 
Jn 

whose time derivative is given by the flux through the 
boundary, 

E~ [ cj) n dS, (76) 
Jan 

where n, is now the outward pointing normal to the 
boundary <9f2. In Appendix [21 we show that symmet- 
ric hyperbolicity implies strong hyperbolicity and that 
the boundary flux can be written as 

0»=5>al4. (77) 

a 

where the sum is over a suitable basis of characteristic 
variables. Therefore the growth of the energy can be con- 
trolled by controlling all characteristic variables that are 
ingoing at the boundary (Aq, > 0) ("maximally dissipa- 
tive boundary conditions"), while ingoing or tangential 
characteristic variables give negative or zero contribu- 
tions to 4> n and hence E. 

The same result holds for a second-order system but, as 
we have seen, the characteristic variables of the second- 
order system are not unique: the v t A can be given arbi- 
trary speeds, and arbitrary combinations of v t A can be 
added to any characteristic variable. In order to impose 
maximally dissipative boundary conditions with the de- 
sired effect of controlling an energy, we need to control all 
incoming characteristic variables of an actual first-order 
reduction. We can do this within the second-order sys- 
tem by imposing boundary conditions on characteristic 
variables u Q constructed from w and v j, but we need to 
have the correct admixtures of v t A i n these characteristic 
variables. 

We have seen that in order to show strong hyperbolic- 
ity for the second-order system one needs to diagonalise 
only the matrix A. Let us call the characteristic vari- 
ables of the second-order system given by eigenvectors of 
A\ and which are therefore constructed only from by w 
and V !n , the "short" characteristic variables, and let us 
denote them by u^, . For proving symmetric hyperbolicity 
one only needs to find a conserved energy e(w, u j). 

Assuming that we already have the u^, on the one 
hand, and the energy e(w,v t i) on the other, the simplest 
way of generating the required "full" characteristic vari- 
ables u Q is to make the ansatz u Q = const (u' Q + un- 
determined multiples of the d^), and to determine the 
coefficients and overall normalisation for each u a so that 
the full characteristic variables obey 

e = 5> 2 a . (78) 

a 

(|76|l then follows. 



V. PSEUDO-DIFFERENTIAL REDUCTION 
METHOD 

For the purpose of comparison, we now describe the 
pseudo-differential reduction method of 0, 13 in our no- 
tation. We carry out a Fourier transform in space with 
wave number u>i of the second-order system <|11I12[1 . We 
denote the Fourier transforms of v and w by v and w. 
We choose the direction ni to be that of the wavenum- 
ber oji and write LJi = \uj\rii. We then have d n = i\uj\v 
and d,A = 0. We can then use d n to represent v. The 
principal part of the Fourier-transformed system can be 
written as 

where A is the matrix defined above, and the non- 
principal terms not shown here are homogeneous of order 
zero in |u;|. In our notation, the definition of § is then 
Definition 2c: The second-order system is called 
strongly hyperbolic if and only if there exists a pseudo- 
differential reduction to first order of the form Tfflti where 
A is diagonalisable with real eigenvalues, and the diago- 
nalising matrix and its inverse depend smoothly on rii. 

We have set up our notation so that from Theorem 1 
we immediately have 

Lemma 5: Definition 2c is equivalent to Definition 2a 
and to Definition 2b. 

As the pseudo-differential approach relies in an essen- 
tial way on Fourier transforms, it does not lend itself to 
defining a locally conserved energy. Therefore there is no 
definition of symmetric hyperbolicity in this approach. 



VI. THE EVOLUTION OF CONSTRAINTS ON 
THE ORIGINAL SECOND-ORDER SYSTEM 

In many applications, the original second-order sys- 
tem is subject to differential constraints, which are con- 
served under evolution. We shall call these the "original 
constraints" to distinguish them from the "auxiliary con- 
straints" di — v t i — that arise additionally in the process 
of first-order reduction. 

Note that until now we have not mentioned or used 
original constraints. The reason is that in general one 
wants to prove well-posedness of the second-order sys- 
tem if the original constraints are obeyed or not. This is 
important for example if one wants to carry out numer- 
ical simulations using "free evolution" where the origi- 
nal constraints are imposed only on the initial data. In 
the continuum the constraints then remain zero, but in 
numerical free evolution they are violated through finite 
differencing error. At a later time one is effectively evolv- 
ing initial data that do not obey the constraints. There- 
fore the continuum problem must be well-posed for non- 
vanishing original constraints as a necessary condition for 
numerical stability. One may of course use the original 
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constraints to modify the original, second-order system, 
but here we assume that this has already been done, and 
the second-order system is fixed. 

We now prove that the evolution system of original 
constraints is strongly hyperbolic if it closes and if the 
second-order main system is strongly hyperbolic. An 
equivalent result for first-order systems subject to con- 
straints is given in |13J. 

We consider a vector c of constraints which are quasi- 
linear of the form 

c ~ Cfv^ + C % 2 w,i, (80) 

where the matrix C lJ is symmetric in ij. It is clear from 
(|11I12[) that the evolution of these constraints is first or- 
der in space and time. If the constraint system is closed, 
its principal part must then be of the form 

c~G l c :l (81) 

for a vector of square matrices G % . Using the second- 
order evolution equations H11I12J) and comparing the 
leading order terms in v and w we find 

C^A^+C ( 2 B{ k) = G^C{ k \ (82) 
C^Ai + CgBg = G (l C ] 2 \ (83) 

These identities between totally symmetric matrices hold 
if and only if their contraction with m on all indices hold 
for all rii- Writing the nnn and nn components of these 
equations in matrix form we have 

/ ™ r>nn \ 

(c ,„ CT) I ^2 ^ \ = G»(c»,cr) (84) 

for all rij. We write this in compact notation as 

CA = GC. (85) 

If the second-order main system is strongly hyperbolic, 
A is diagonalisable with A = TAT -1 . G can always be 
brought into Jordan form as G = SJS -1 . Then 

CA = JC, C = S- l GT, (86) 

We assume that the rows of C, and therefore the rows 
of C, are linearly independent. This means that there is 
no redundancy between the differential constraints, and 
is similar to an assumption in [l3| . 

Consider now the first Jordan block of J with eigen- 
value \x\. For simplicity assume it has size 2. Exception- 
ally writing out the internal matrix indices, we have 

G\a\ a — l^lGia + C 2a , 

C^aXa = f-lC 2a (87) 

(no sum over the index a). From the second equation 
C*2a = for all a such that \ a ^ Using this result 
and the first equation, C 2a = precisely for those a for 
which X a — fii. By assumption, no row of C vanishes, 



so C 2a cannot all vanish. Therefore, there must be at 
least one a such that X a = /j,x, and the first equation 
must be absent, that is, the Jordan block has only size 
1. Repeating this argument for all Jordan blocks of J 
means that each eigenvalue of G coincides with one of 
A, and that J is diagonal, that is, G is diagonalisable. 
Writing G — SA'S^ 1 where the diagonal matrix A' is a 
submatrix of A, we have 

(SC)A = A'(SC). (88) 

This means that the rows of SC are left eigenvectors 
of A, and parameterise to characteristic variables of the 
second-order system. We have shown 
Theorem 5: The evolution of the original constraints 
is strongly hyperbolic if the second-order main system is, 
and its characteristic speeds are then a subset of those 
of the main system. Furthermore, we can find a basis 
of characteristic variables for the main system and the 
constraint system such that for each characteristic vari- 
able c a of the constraint system, there is a characteristic 
variable u a of the main system such that 

c a = d n u a + transverse derivatives. (89) 

Note that there is no such result for symmetric hyper- 
bolicity. 

VII. SYMMETRIC HYPERBOLIC-PARABOLIC 
SYSTEMS 

Theorem 4.6.2 of asserts the following: Assume we 
have a vector of variables u and another vector of vari- 
ables z, which obey a linear system of evolution equations 
of the form 

d t u = D n u + D 12 z, (90) 
d t z = D 21 u + D 22 z. (91) 

Here the D are linear spatial derivative operators whose 
coefficients can depend on t and x l . Dn is a first-order 
derivative operator such that dtu — D\\u is symmet- 
ric hyperbolic. D 22 is a second-order derivative operator 
such that dtz = D 22 z is parabolic. D\ 2 and TJ21 are arbi- 
trary first-order derivative operators. Then the coupled 
system is called mixed symmetric hyperbolic/parabolic. 
Its Cauchy problem with periodic boundaries is well- 
posed. 

The theorem can be applied straightforwardly to 
second-order systems. We identify the variables u of the 
theorem with the variables (v, w, d{) of the first-order re- 
duction of what is to be the symmetric hyperbolic sub- 
system, and then go back to the second-order form of 
this subsystem by replacing di with . The result is the 
system 

v = A\vj. + A-, v + A->.w + a + Cz, (92) 
lb = B\ 3 v^j + B[v^ + B x v 4- B\w t i + B 2 w + b 
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+D l z. l + Dz, (93) 

z = D 2 2Z 

+ Ej j v tij + Ejv lt + E lV + Ejw,j + E 2 w. (94) 

The coupling operators D12 and D21 are parameterised 
by the matrices C and D, and E, respectively. We have 
underlined their principal parts to show what order of 
derivative is allowed in the coupling terms. 
Definition 6: A second- order system is called mixed 
symmetric hyperbolic-parabolic if it is of the form 
\y$ , such that D22 is parabolic and the system ^HHlfy) 
with the same coefficients is symmetric hyperbolic (in the 
sense of Definition 2a or 2b). 

Theorem 4.6.2 of [l| then gives us 
Lemma 6: The Cauchy problem with periodic boundary 
conditions for such a system is well-posed. 

VIII. CONCLUSIONS 

We have formalised the definition of strong or symmet- 
ric hyperbolicity of a system of evolution equations that 
are first order in time and second order in space by re- 
ducing them to an equivalent first-order system. We have 
given necessary and sufficient criteria for the existence of 
a reduction that is strongly hyperbolic or symmetric hy- 
perbolic. These criteria are formulated entirely in terms 
of the principal part of the second-order system, without 
an explicit reference to the reduction. 

We have proved that the definitions of strong hy- 
perbolicity based on a first-order reduction, a pseudo- 
differential reduction, and a direct second-order approach 
are all equivalent. The definitions of symmetric hyper- 
bolicity based on a first-order reduction and on a direct 
second-order approach are also equivalent. 

In order to analyse the well-posedness of a given 
second-order system in practice, there are three non- 
trivial calculations to be carried out, independently of the 
approach in which one has defined hyperbolicity. Sup- 
pressing technical details, they are as follows. 

Strong hyperbolicity Strong hyperbolicity of the 
second-order system is equivalent to diagonalisability, 
with real eigenvalues, of the matrix A. 

Symmetric hyperbolicity Symmetric hyperbolicity of 
the second-order system is equivalent to the existence of 
an energy and flux quadratic in the and w. 

Maximally dissipative boundary conditions In order 
to impose maximally dissipative boundary conditions, 
one needs the full characteristic variables of a symmetric 
hyperbolic first-order reduction. This is done most easily 
starting from the left eigenvectors of A and the energy 
e(w,v ti ). 

We have established criteria for well-posedness of the 
second-order system regardless of any constraints it is 
subject to. However, if the second-order system is 
strongly hyperbolic, and there is a closed constraint sys- 
tem associated with it, then the constraint system is also 
strongly hyperbolic, and the characteristic variables of 



the constraint system are related to a subset of the char- 
acteristic variables of the main system. 

It is known that a first-order symmetric hyperbolic sys- 
tem coupled to a parabolic system through at most first 
derivatives of all variables has a well-posed Cauchy prob- 
lem. We have generalised this result to second-order sym- 
metric hyperbolic systems through a reduction to first 
order. 

Appendix [D] gives an example of a simple second or- 
der system discussed in both the first order reduction 
approach and the direct second-order approach. 
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APPENDIX A: DIAGONALISABILITY OF 
MATRICES WITH A ZERO BLOCK 

Consider the matrix 

where the square block A has size n and the square block 
C has size m. We now prove that if M is diagonalisable 
then A and C are both diagonalisable. 

The eigenvalues of M are clearly the union of those 
of A and C . The eigenvectors of M can be constructed 
from those of A and C as follows: Suppose we have 

Avi = XiVi, i = l,...,n, (A2) 
Cuij = ^jWj, j = l,...,m. (A3) 

Then we can form the eigenvectors Xi — (vi,0), which 
span the invariant subspace of M, and yj = (uj, Wj) such 
that 

Mx l = XiXi, i = l,...,n, (A4) 
M V] = NVh = 1, • • • , m, (A5) 

with the condition (A — [ij)uj = —Bwj for j = 1, . . . , m. 
If /ij is not an eigenvalue of A then (A — /ij ) can be in- 
verted and there is a unique solution for Uj . Therefore, if 
the eigenvalues of A and C are disjoint, the eigenvectors 
can be completed and the linearly independent 

from the yj because they correspond to different eigen- 
values. In this case M is diagonalisable if and only if 
both A and C are diagonalisable. 
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Now suppose that A and C share an eigenvalue A. The 
transformation 



From this and 



(A6) 



with 



U = 



Un 







(A7) 

and D = UxiXC-AX + Bpz 1 brings 4 and C into Jor- 
dan form simultaneously for suitable U\ and 11% . Without 
loss of generality we can assume that A is a single Jor- 
dan block of eigenvalue A and rank(A — XI) = r, and 
that C is another Jordan block with the same eigenvalue 
and rank(C — XI) = s. The matrix M is diagonalisable 
if rank(Af - XI) = 0. We have rank(M - XI) > r + s 
because the r + s columns of (M — XI) containing a 1 in 
the second diagonal are linearly independent, while the 
matrix B could provide additional linearly independent 
vectors. Therefore if M is diagonalisable we must have 
r = s = for each Jordan block of A and C, and so A 
and C are diagonalisable. 



APPENDIX B: HYPERBOLICITY OF THE 
AUXILIARY CONSTRAINT SYSTEM 

With the choice of reduction parameters that we have 
used in the proof of Theorem 1, namely l|29|) . (|30[1 . and 
using the following further auxiliary constraints, 



Cij Cij -\- C[z,j] 0, 



Cijk 



Cii.k i Cjk.i i Chi, 



o. 



(Bl) 
(B2) 



the auxiliary constraint system can be reduced to the 
decoupled system 



Aid - i^ei jk Cj^, 



A x c i3 - ine[i kl c k i t \j] 



(B3) 
(B4) 



This is strongly hyperbolic with speeds 0, ±/i for all 
Furthermore 



Ci = AiCi-inu^C^k, 



(B5) 
(B6) 



which is also strongly hyperbolic (in fact, symmetric hy- 
perbolic), and so in this case the auxiliary constraint sys- 
tem is well-posed. 



APPENDIX C: SYMMETRIC HYPERBOLICITY 
AND CHARACTERISTIC VARIABLES 

Assume that P % admits a symmetriser H . As if is Hcr- 
mitian and positive definite, there is an invertible matrix 
S such that 



H = S^S. 



(CI) 



HP' = P l] H 



(C2) 

it follows that SP n S~ 1 is Hermitian, for any direction 
rii. Therefore it can be diagonalised by an orthogonal 
matrix R (which generally depends on rii), or 



SP n S~ 1 = R-^AR, 
where A is diagonal. Therefore 

pn = T -l ATj rp = 



(C3) 



(C4) 



and we have proved that symmetric hyperbolicity implies 
strong hyperbolicity. 

Furthermore, as R is orthogonal, 



H = S\R)R)S — T^T, 



(C5) 



and so there is a preferred basis, namely the rows of T 
(which generally depends on rii), of left eigenvectors of 
P n in which H is the unit matrix. In terms of the original 
basis 

H = T^T, HP n = T f AT. (C6) 
In quadratic forms the same fact can be written as 

a a 

where the sum is over the characteristic variables in the 
basis encoded in the rows of T. 



APPENDIX D: THE KWB FORMULATION OF 
THE MAXWELL EQUATIONS 

We use this formulation of the Maxwell equations to 
illustrate some of the differences between the reduction 
approach and the second-order approach, namely the ex- 
istence of "short" and "full" characteristic variables, the 
split of the free parameters of H into those that are re- 
duction parameters and those that are not, and the rela- 
tion between the characteristic variables of the main and 
constraint systems. We work initially in the first-order 
reduction approach, and then re-interpret the same cal- 
culations in the language of the second-order approach 
afterwards. 

The system has been discussed in 0, 0|. In flat 
spacetime, in radiation gauge and in the absence of 
charges (source terms), it is 



Ai — E{ , 



E, 



-Aij' j + (1 - -/i.l , , ' • <il... 



(Dl) 
(D2) 
(D3) 

where v = {Aj}, w = {Ei,T} are the dynamical vari- 
ables, and a and b are constants that parameterise addi- 
tion of the ( "original" ) constraints 



r = (b-i)E\i 



Ce 



r - as = o, 

E/ = 



(D4) 
(D5) 
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to the evolution equations. Repeated indices have been 
raised with the metric S lJ and are summed over. We 
shall also use 8ij to decompose tensors into their normal 
and transverse parts. With dij = Ajj, we consider the 
first-order reduction. 

A t = -E u (D6) 
Ei = -dji' j +aT ti + (l-a)d j j i 

•'('/./' '/-'.,)• (D7) 
f = (b-l)E\ h (D8) 
dij = -E jti . (D9) 

The constant c parameterises , and ZV fe has been 
set to zero. For the other reduction parameters we have 
made the standard choice H29|) . so that the Ai decouple 
from the Ei and T. 

The matrix V , leaving out the zero rows and columns 
corresponding to Ai , is block-diagonal with the blocks 



a —a 1 — a — c 
b-1 
-10 





(o) (d AB y (D12) 



Here d qq = q lJ dij, and d A B represents the transverse 
trace-free object qfq l jdki — (l/2)qijd qq . The characteristic 
variables are 

u = T + (b-l)d nn , (D13) 
u± = a(T-d nn ) + (l-a-c)d qq ±VabE n (pU) 
u± A = d nA - cd An T E A , (D15) 



with speeds (0, ±vafe, ±1), and d qq , d An and d A B with 
zero speed. 

The first-order reduction admits the conserved energy 



e = E i E i + d. l0 d l] - 2aTd i i + (2a - 1 - ab)(d 



» ^2 



+ci [r + (6-iK i ] + c [(d*y-d«(p'] (Die) 

with the flux 



2 [o(r - d j r )E i + ,1, ' /:' - <f J '£j-] 
+2c(d? i E j -dj i E i ), 



(D17) 



where ci is a free parameter in the energy, and c is the 
reduction parameter introduced above. 

We now review how one would deal with the same sys- 
tem in a direct second-order approach [pointing out the 
relation to the first-order approach in square brackets] . In 



order to show strong hyperbolicity, one would diagonalise 
the matrix A. It is block-diagonal with the blocks 



(D18) 




and 



-1 



E A 

A A , n 



(D19) 



[A is the sub-matrix of V obtained by suppressing the 
rows and columns relating to d A {. In a different point 
of view, we could set the rows and columns relating to 
d A i in V to zero by allowing the reduction parameters to 
depend explicitly on In our example, this corresponds 
to setting c = 1 - a in (|D10fl but c = in (jDll|) .] 

The "short" characteristic variables of the second- 
order system are obtained as eigenvectors of A'. They 
are 



U n 



= r + (6-i)A„, n , 



u' ± = a(r — A n<n ) ± VabE„ 



l ±A 



A A , n T E A , 



(D20) 
(D21) 
(D22) 



with speeds (0, ±\/a&, ±1)- [These are the characteristic 
variables of the reduction, minus all terms in d A i = Ai tA .] 
To show symmetric hyperbolicity of the second-order 
system, we find an an energy. The most general one is 



e = E,E l 



2aTAi l + (2a - 1 - ab)(A l ' 1 ) 2 



+ci [T+(b-l)A/] +c 2 [(A t 



A^jA" 1 ' 



(D23) 



with the flux 



2 [a(T 
+2c 2 (A i ' j E j 



A 3 ' 3 )E l 



A 



i r A j ' j E i - A^Ej] 
' j E l ), (D24) 



where c\ and C2 are free parameters. [This is identical 
to the energy of the reduction, except that c 2 is now not 
a reduction parameter. Rather, the term it multiplies 
is independently conserved if we allow commutation of 
partial derivatives.] 

In order to impose maximally dissipative boundary 
conditions, one needs the full characteristic variables u Q . 
We find this by expressing e as a quadratic form in u Q . 
Comparing e with the suggests that 



1 

2ab 



+ i (u +J 4U +yl + u_ A u A ) 

n - 7 j U5 - c 2 (u +A + u- A )A n ' A 
quadratic in A i A , (D25) 



(«■-!) 



with u a = u' a + multiples of the A i A . [The result is 
equivalent to the u Q of the reduction, with dij — ► Aj^ 
and c — > C2-] 
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Finally, the constraint evolution system is and these expressions hold for both the second-order sys- 

tem and the reduction. 

C r = bC E , (D26) 
C E = aC T /, (D27) 

Its non-trivial characteristic variables are 

c± = ad„Cr ± VoECe = <9„u± + transv. deriv., (D28) 
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